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Abstract 

Using path-integral Quantum Monte Carlo we study the low-temperature 
phase diagram of a two-dimensional superconductor within a phenomeno- 
logical model, where vortices have a finite mass and move in a dissipative 
environment modeled by a Caldeira-Leggett term. The quantum vortex liq- 
uid at high magnetic fields exhibits superfluidity and thus corresponds to a 
superinsulating phase which is characterized by a nonlinear voltage-current 
law for an infinite system in the absence of pinning. This superinsulating 
phase is shifted to higher magnetic fields in the presence of dissipation. 
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Quantum fluctuations of vortices play an important role for the low-temperature physics 
of two-dimensional (2D) superconductors. It has been established that dissipative Jy], inertial 
PJ3[], or Hall-type M zero point motion is able to melt the vortex lattice at sufficiently high 
magnetic fields and give way to a quantum vortex liquid phase in principle at T = 0. In 
the presence of disorder this melting transition is a continuous one and corresponds to the 
localization or crystallization of Cooper pairs in the theory of the superconductor-insulator 
transition However, for a system without disorder, the melting of the vortex lattice is 

discontinuous in (2+1) dimensions JJ. 

In this letter, using path-integral Quantum Monte Carlo (QMC), we will study a phe- 
nomenological model where the vortices have a finite mass and move in a dissipative envi- 
ronment, which is modeled by a Caldeira-Leggett term. Both, vortex mass and dissipation 
are assumed to originate from electronic contributions while the vortices themselves are con- 
sidered as bosons. For the simulations we focus on the situation where the Magnus force 
is zero which is appropriate for the case of granular systems or Josephson Junction arrays 
(JJA) where experiments suggest that the transverse force is very small The effect of a 
weak Magnus force will then be discussed afterwards. 

In a thin film superconductor at high magnetic fields vortices effectively interact loga- 
rithmically. So the question arises whether the vortices can form a superfluid or not. We will 
show using QMC that a 2D system of logarithmically interacting particles indeed exhibits 
superfluidity at low temperatures and high densities even in the presence of dissipation. || 
This means that the quantum vortex liquid actually corresponds to a superinsulating phase 
in which an infinitesimal current will cause an infinite voltage in the absence of pinning effects 
in analogy with sliding charge density waves. This superinsulating phase exhibits a non- 
linear volt age- current law with non-universal exponent and is different from the insulating 
Cooper pair glass phase in Fisher's || treatment of the field-tuned superconductor-insulator 
transition in the presence of strong disorder. We will show numerically that at finite tem- 
perature the superinsulator gives way to a classical vortex liquid via a Kosterlitz-Thouless 
(KT-) transition. In the presence of dissipation this phase transition is shifted to higher 
magnetic fields possibly giving rise to an intermediate non-superfluid vortex liquid phase at 
T = 0. 

In the non-dissipative case the vortex system is described by the Euclidean Lagrangian 

where rfe(r) are the vortex positions, m is the vortex mass, and To controls the strength 
of the logarithmic interactions. We assume the system to be embedded in a neutralizing 
background which corresponds to the presence of a uniform external magnetic field B = $ f? 
so that the long-range interaction in Eq. (H) is well-defined even for periodic boundary 
conditions. Here, g is the vortex density, and $o is the flux quantum. r^(r) obeys periodic 
boundary conditions in imaginary time r, that include bosonic exchange of the particles. 

Finite-temperature QMC simulations have been performed for systems with up to 128 
time slices and N=16, 28 and 36 particles on a grid of 256x222 sites in a rectangular periodic 
box which was commensurable with a triangular vortex lattice. A bisection algorithm was 
used where several particles are updated in several imaginary time slices simultaneously in 
order to allow for cutting and reconnection of the paths [Q. The superfluid density g s was 
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then obtained from the distribution of winding numbers W;, (i = x,y) around the periodic 
cell by @ 

Q s rnL] 



Wf) , (2) 



g (3h 2 N 

where Li is the periodic box length in i-direction, and /3 is the inverse temperature. For the 
lowest-temperature data reported in this work typically about 10 5 sweeps were needed to 
equilibrate the winding number moves. Equilibration was checked by carefully monitoring 
the distribution of Wi which can be fitted by a Gaussian at low T. In the following we use 
dimensionless variables, where the temperature T is given in units of To, and the magnetic 
field B is given in units of B = (mT Q )/(y/3 h 2 ). 

The phase diagram of the model (HD is shown in Fig. la. It consists of three phases: a 
superconducting vortex solid phase (VS) at low densities and temperatures, a classical vor- 
tex liquid (VL) at high temperatures, and a quantum vortex liquid (QVL) at high densities 
and low temperatures. The vortex solid thermally melts via a KT-transition mediated by 



the unbinding of dislocation pairs which takes place at Tkt ~ 0.0071 |16j and is approx- 
imately independent of the magnetic field B. At low temperatures quantum fluctuations 
melt the vortex solid via a discontinuous transition at the melting field Bm ~ 0.0077 
which is approximately independent of the temperature . In order to check our numerical 
algorithm we measured the internal energy at various fields and low temperature, and also 
estimated Bm by calculating the Bragg-peak intensity. The results where found to be in 
good quantitative agreement with the zero temperature results reported in Ref. [0]. 

The inset in Fig. la shows the superfluid fraction g s / g as function of the magnetic field 
B at T = 0.005 as obtained in our simulations. At the quantum melting transition g s /g 
jumps from zero to one which quantitatively reproduces the result of Nordborg and Blatter 
H, and shows that the QVL actually is a vortex superfluid with g s = g. In the following 
we will focus on the transition from the QVL phase to the VL phase. The behavior of the 
superfluid fraction is shown in Fig. 2a where g s / g is plotted as function of the temperature 
T for different magnetic fields B. It is seen that all data collapse onto one single curve 
when plotted against the scaling variable T/B. This scaling behavior can be understood 
with a KT-type unbinding of topological (vortex-antivortex) excitations (which we will call 
dual vortices) in the vortex-superfluid. Since the energy scale of the logarithmic interaction 



of dual vortices is given by Tq = 2nh g/m, the transition temperature thus scales with B. 
Because of the fundamental duality between vortices and charges |K| the dual vortices are 
actually related to the Cooper pair degrees of freedom in the system as will be discussed 
below. 

Following Ceperley and Pollock [|l7j] further numerical evidence for the existence of a 
KT-transition can be obtained by performing a finite-size analysis that explicitly invokes 
the KT-recursion relations |18| which are integrated up to the system radius L/2. There 
are only two independent fit parameters for all numerical data involved in this procedure: 
the dual vortex core energy E and the dual vortex radius d. The resulting fits are the solid 
curves in Fig.s 2a and 2b. Fig. 2b shows data for iV =16, 28, 36, and also extrapolated 
to iV = oo. From this the KT-transition temperature is obtained as T' KT = 1.451 h 2 g/m, 
where E = 6.49 h 2 g/m, and d = 0.634 g' 1 / 2 . Interestingly it is observed that the core size 
d scales with the vortex spacing g~ x l 2 which is a consequence of the logarithmic interactions 
which provide no additional length scale. 
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In order to interpret these numerical results we note that the vortex superfluid actually 
should be considered as a charged superfluid where the flux quantum $ plays the role of 
a "charge" and the "flux quantum" of the dual vortices is given by 2e. In a homogeneous 
superconductor the corresponding gauge field a is related to the 2D-superfluid density of 
Cooper pairs n s , V x a = 2en s , which gives rise to the Magnus force acting on the vortices 
||1 1|| . To be consistent with the assumption of of a zero transverse force in a granular 
system we therefore suppose that n s is replaced by the density of fluctuating charges 5n s 
on the superconducting islands, V x a = 2e5n s , which is a usual assumption in the context 
of granular superconductors or JJA's [T2|]. Static fluctuations of the gauge field a can 



approximately be included by a contribution (A/4e 2 )(V x a) 2 to the free energy of the 
vortex superfluid, where A is the strength of the (S-) interaction between the Cooper pairs 



22]] . The numerical simulations then correspond to the case of strong repulsion, A = oo, 
where 5n s decays logarithmically around a dual vortex. For finite A the vortex-superfluid 

current flowing around the dual vortices is screened on the length A = (2Am/h 2 g^J ^ . In 
the following we assume that this screening length is much larger than the dual vortex core 
size d. 

Let us discuss the dynamics of the vortex superfluid in the case A = oo. Since dual 
vortex excitations carry a charge ±2e the separation of dual vortex-antivortex pairs leads 
to an electric current J. However, in an infinite system this process requires an infinite 
amount of energy when the electric field E is zero, while for a non-zero electric field the 
energy barrier is finite. It is well known [TJ3| that for T <C T' KT this type of dynamics leads 



to a non-linear voltage-current law with nonuniversal exponent, 

E~J T / T o. (3) 

Eq. explicitly shows that the vortex superfluid phase actually is a superinsulating phase 
with infinite resistivity even at non-zero temperature. In a realistic system finite-size effects 
as well as a finite screening length A will restore a finite resistivity p ~ e E o/ T for small 
currents where the Arrhenius factor e~ E °^ T controls the number of dual vortex-antivortex 
pairs in the system, and the KT-transition will be broadened. 

Let us now briefly consider the case where a weak Magnus force is present. This means, 
that the vortices on average pick up a phase <fi = 2im° s ^Q 7^ when encircling a supercon- 
ducting grain of area £q |I3| . Here, 2en° s is the charge density of Cooper pairs contributing 



to the transverse force which plays the role of a fictitious external magnetic field. By anal- 
ogy with an ordinary superconductor we therefore expect that dual vortices will form an 
Abrikosov lattice. Thermal excitations and dynamics of this dual vortex lattice will then 
be governed by the unbinding of dislocation pairs and the KT-melting temperature is given 
by T'kj, ^0.06 h 2 g/m. This again is a superinsulating phase described by the non-linear 
volt age- current law (^), where Tq is replaced by the interaction energy of dislocations. The 
corresponding phase diagram is schematically shown in Fig. lb. 

In what follows we will discuss the influence of dissipation, which is treated approximately 
by adding a Caldeira-Leggett term (time-delay) |20| 



hf3 



CcL = -2^Vl dT '^ ln 



sin ^(r-r 



— (4) 
dr> V 1 
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to the Lagrangian (|T|). Eq. phenomenologically describes the coupling of each single 
vortex to a separate harmonic oscillator heat bath which is characterized by a frictional con- 
stant 77. This model is appropriate to describe dissipation originating from normal currents 
in the vortex cores [[J and has also been used by Wheatley |2l| to study the statistics of 
holons coupled to the background of spinons in the context of high-T c superconductivity. 
With regard to superfluidity there is one essential assumption in using Eq. to model a 
dissipative environment. The coupling to independent harmonic oscillator heat baths can- 
not exactly be valid for indistinguishable particles. Following we therefore assume that 
dependencies between separate heat baths induced by bosonic exchange of the vortices can 
be neglected. 

We again computed the superfluid fraction which is shown in Fig. 3 for B = 0.04, 0.08, 
and 7] = 0.01, 0.02, where r] is given in units 7/0 = TtmTo/h. It is seen that dissipation 
leads to a shift of the KT-transition to lower temperatures and, more surprisingly, to the 
suppression of superfluidity at T = 0. This reentrance behavior can be understood by 
the contribution of the Caldeira-Leggett term to the action of vortices involved in a 
particle exchange, which can be described in terms of a temperature-dependent effective 
mass, m* = m + /3hri(\ia2 j V) . For the bare vortex superfluid density g s (unrenormalized by 
dual vortex-antivortex pairs) we therefore make the ansatz g s = g(m/m*), which is a linear 
function of T at low temperatures. Since the interaction of dual vortices scales with g s , the 
KT-transition is then shifted to lower temperatures, Tj^ T = T' KT — (rjh / 'm) (In 2 /tt). Using 
this ansatz and taking the core-energy to be Eq = E$q s / q we perform the same finite-size 
fit as above which yields the solid curves in Fig. 3 and shows a remarkable coincidence with 
the numerical data. Again the same values for Eq and d as before were used so that no 
new fit parameters appear here. The dashed curves in Fig. 3 show the superfluid fraction 
extrapolated for N = 00. At low temperature, g s /g rises linearly with T and reaches its 
maximum value ~ 1 — 0.15 r//(hg) near T'^ T . 

The KT-transition line between the VL and the QVL (extrapolated to T = 0) is shown 
in the S-T-diagram in the inset of Fig. 3 where the friction parameter rj was assumed to 
be independent of the temperature. The crossing point with the I?-axis (T^t = 0) * s given 
by B Q KT « 0.15 ri& /h. If B° KT > B M this suggests the existence of an intermediate field 



range in which the classical vortex liquid phase continues to T = 0. |23[ In this regime the 
superinsulating phase is unstable against the unbinding of dual vortex-antivortex pairs even 
at T = since the strength of the logarithmic dual vortex-vortex interaction is proportional 
to T. This also means that the exponent of the voltage-current law will stay finite in 
the limit T — > 0. 

We add a few comments about the experimental accessibility of the results presented in 
this paper. The notion of a quantum melting transition as it was treated here requires a 
system with weak pinning as well as the existence of a finite vortex mass which has to be 
small enough in order to push the melting field Bm significantly below the upper critical field 
B c2 . Indeed, a field-tuned superconductor- insulator transition which can be related to the 
melting of the vortex lattice has been observed in Josephson junction arrays [23]. According 
to our results for a finite system the conductivity in the superinsulating phase is expected 
to show an Arrhenius behavior a ~ e - E o/T w h ere dissipation (as it was modeled here) leads 
to a linear T-dependence of E at low temperatures so that the conductivity stays finite in 
the limit T — > 0. This is qualitatively consistent with the experimental results in Ref. pi. 
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In conclusion, using QMC simulations, we have studied superfluidity in a system of log- 
arithmically interacting massive vortices moving in a dissipative environment. It was shown 
that the quantum vortex liquid at high magnetic fields corresponds to a superinsulating 
phase characterized by a non-linear voltage-current law with non-universal exponent. This 
superinsulating phase gives way to a classical vortex-liquid phase via a Kosterlitz-Thouless 
transition, which is shifted to higher magnetic fields in the presence of dissipation. 
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FIGURES 




FIG. 1. (B-T)-phase diagram of the vortex system in a 2D superconductor (below B C 2). The 
inset shows the jump in the superfluid density at the quantum melting transition at % ps 0.0077 
(arrow), (a) zero Magnus force, (b) weak Magnus force. 
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FIG. 2. Superfluid fraction q s /q: (a) as function of the scaling variable T/B for the magnetic 
fields B = 0.03, 0.04, 0.08, and 0.12, and a 28-particle system, (b) at B = 0.04 as function of T 
for systems with 16, 28, and 36 particles and extrapolated for an infinite system. 
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FIG. 3. Superfluid fraction g s /g as function of T for non-zero dissipation (N = 28): B = 0.04, 
r) = 0.01 (O), B = 0.04, n = 0.02 (A), and B = 0.08, rj = 0.01 (□). Dashed curves: Finite-size 
extrapolation to N = oo. Inset: (S-T)-phase diagram (schematically). 
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